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I review domain wall fermions in vector gauge theories. Following a brief introduction, the status of lattice 
calculations using domain wall fermions is presented. I focus on results from QCD, including the light quark 
masses and spectrum, weak matrix elements, the n/ = 2 finite temperature phase transition, and topology and 
zero modes and conclude with topics for future study. 



1. INTRODUCTION 

Several years ago Kaplan had the idea for a lat- 
tice discretization of the Dirac operator that pre- 
serves chiral symmetry at any lattice spacing [Q. 
This miracle is performed by adding an odd in- 
finite dimension to the usual even dimensional 
spacetime. The mass of the fermion is given the 
shape of a domain wall in this extra dimension. 
The mass is large, positive on one side of the de- 
fect, negative on the other. Chiral zero modes 
then arise naturally on the wall where the mass 
is zero. If the extra dimension is periodic, an 
anti-domain wall appears at the other end of the 
lattice with a zero mode of the opposite chirality. 
These two chiral zero modes are coupled to the 
same (4d) gauge field to simulate a vector gauge 
theory, like QCD. The chiral symmetry is mani- 
fest: the left and right handed fermions have been 
separated globally in the extra dimension, so they 
may be rotated independently. 

This realization of chiral zero modes was al- 
ready known to occur in the continuum, for in- 
stance, by Callan and Harvey who studied topo- 
logical defects embedded in higher dimensions!^. 
There and in Kaplan's original study one may 
suppose that the extra 5th dimension is physi- 
cal: our 4d spacetime is bound to the wall at low 
energies, but at high enough energies the 5th di- 
mension becomes perceivable. It is interesting to 
note such ideas have recently been proposed as a 
natural mechanism for the spontaneous breaking 
of supersymmetryj^. Loosely speaking, domain 
walls are D-branes, and similar models for real- 



istic "compactified" universes have recently been 
proposed(for example, see Ref. Q]). 

Shortly after Kaplan's original proposal and 
also motivated by an independent proposal by 
Frolov and SlavnovQ, Naryanan and Neuberger 
transformed the extra dimension into an infinite 
flavor space of heavy regulator fields Using 
a transfer matrix formalism, the chiral determi- 
nant is given as the overlap of the ground states 
of two Hamiltonians (hence the name, the overlap 
formulation, or the overlap, for short). The over- 
lap and Kaplan's domain wall fermions (DWF) are 
equivalent in the particular case where the extra 
dimension is infinite. The overlap provides an el- 
egant framework for studying the Index theorem 
on the lattice. It is the starting point for a new 
lattice Dirac operator that also maintains chiral 
symmetry on the lattice Q (see below). An exten- 
sive review of the overlap is given in Ref. 

The original proposal by Kaplan was intended 
as a lattice chiral gauge theory. However, it was 
soon recognized that the discretization would be 
useful for lattice QCDj|] since conventional meth- 
ods explicitly break chiral symmetry. Recent nu- 
merical studies have shown the actual extent of 
the extra dimension can be modest, 10-20 sites for 
simulations of quenched QCD, and still preserve 
the chiral symmetry to a high degree of accu- 
racy ||l^-|l^. Further, DWF are accurate to order 
O(a^) since no (mass) dimension five operators 
that are also chirally symmetric exist to cancel 
0(a) errors! 13 1^. Indeed, the initial numerical 
results indicate improved scaling. This scaling (if 
it holds up after further scrutiny) offsets the cost 
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of the extra dimension. 

For a finite extra dimension with Ng sites, the 
leading discretization errors are expected to be 
0{a) exp(— C/Vg) where C is the rate of exponen- 
tial fall off in the extra dimension. This has been 
verified at one loop in perturbation theory for 
QCdQ. In particular, the authors of Ref. 
calculate the one loop fermion self-energy and 
find that the mass of the light mode remains 
exponentially small. Even after accounting for 
exponentially small terms in the free propagator 
that were neg lected in Ref. Q, Ref. also con- 
cludes that the exponential suppression stays in- 
tact after one loop radiative corrections. 

At this meeting impressive new results were 
presented for lattice QCD using DWF, including 
the first dynamical simulations After a brief 
description of the method, I focus on the lattice 
QCD results. 

2. BOUNDARY FERMION VARIANT[§ 

Shortly after Kaplan's discovery, Shamir refor- 
mulated DWF using a lattice in the extra dimen- 
sion that is half as big. The boundaries of the 
5th dimension are explicitly coupled by a small 
parameter, —m{m > 0). We will see that m is 
proportional to the Ad quark mass and receives 
only multiplicative renormalizations in the limit 
A^s ^ oo. The chiral hmit of DWF is thus 
A^s — > oo and then m — )■ 0. These features 
make Shamir's variant of DWF useful for numer- 
ical simulations. 

The action for DWF is essentially an ordinary 
5d Wilson fermion action: 

Sdwf = - ^ {Px,ySs,s' +Ps,s'5x,y)tp, 
x,y.s,s' 

+ (1 - l^.)UU^.(,,y ) + (M - 4)(5,,„ 
Pr5i,s' - mPLSN,-i,s' - (5o,s' 

(1) 

where M is the full 5d mass parameter, s and 
s' denote the extra dimension, and Pr^l = (1 i 



75)72. There are important differences with the 
ordinary Wilson action, however. The links in 
the s direction are set to unity, the relative sign 
between the Wilson term and the 5d mass term 
is opposite to the usual convention, and, as men- 
tioned above, the boundaries in the s direction are 
coupled with a strength m. The boundary con- 
ditions are anti-periodic since the Kaplan model 
is periodic over 2A^s. Finally, the extra hopping 
terms have 75 's in place of 7^'s; there is no chi- 
rality operator in odd dimensions. 

In the free theory for A^^ — 00 and p^, to — > 
it is straightforward to show that the mass of the 
light quark is iriq = mAI (2 — M) and that the 
singular part of the propagator (on one wall) is 

GRis,.s',p) = p^ M(2^M) 

A similar result holds for the other wall. Aside 
from the extra factor M(2 - M)(l - ^1)"+"' , this 
is same form found in the continuum. The free 
propagator is exponentially suppressed in the ex- 
tra dimension with a rate that depends on M. 
For M = 1 the mode does not penetrate the ex- 
tra dimension at all. Note, the doublers have 
been removed for the choice < AI < 2. As 
M is increased in increments of two, up to ten, 
the zero modes disappear from the spectrum and 
new ones appear, depending on which corners 
of the Brillouin zone, according to the factor 
b{p) = 1 — M + ~ cos(p^)), contribute(see 

^ for details). After M = 10, no zero modes ex- 
ist. The critical value where the light quark first 
appears is Mc = 0. In Ref. the leading order 
corrections to the quark mass for finite Ng were 
given, 

TO, = M{2 - M){m +{1~M + 0(p2))^=). (3) 

The overlap of the two modes in the extra dimen- 
sion has generated an exponentially small quark 
mass. This mass can be made arbitrarily small 
compared to to by increasing Ng. 

In the interacting theory M is additively renor- 
malized, just like ordinary Wilson fermions, M —^ 
M = M — Mc- At one loop, the dominant effect 
comes from the tadpole diagram which is diago- 
nal in hence the above free field results are 
unchanged up to a simple shift in AI . If one uses a 
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reasonable choice for the couphng constant, then 
the magnitude of the tadpole contribution agrees 
nicely with the original nonperturbative estimate 
(near quenched (3 w 6.0) for the optimal value of 
Af which minimizes the overlap of the two light 
modes |]l0| , pl| . A similar shift occurs for the dy- 
namical simulations, as well||l^. Eq. also 
suggests why we do not expect the usual "excep- 
tional" configurations for DWF for even Ng'. the 
contribution to niq is always positive. 

Four dimensional quark fields are constructed 
from the five dimensional fields by taking their 
chiral projections on the boundaries. 



= ^...-i^^+^.o^^^. (4) 

These definitions are the simplest choice for in- 
terpolating operators that create and destroy 
quarks. Other definitions are possible; i.e. one 
may average over some width around each bound- 
ary. 

Operators constructed from the definitions in 
Eq. (I) satisfy a set of axial Ward identities [jisll. 

A,{A;ix)0{y)) = 2m{J^{x)0{y)) 

+i{S%Oiy)) + 2{j^q{x)0{y)). (5) 

The identities are derived in the usual way, by de- 
manding invariance of expectation values under 
infinitesimal transformations. Here, we take ad- 
vantage of the fact that the left and right handed 
modes are globally separated in the fifth dimen- 
sion in order to rotate them independently. The 
first two terms on the r.h.s appear in the con- 
tinuum while the last term is "anomalous" . For 
the fiavor nonsinglet case it vanishes in the limit 
Ng — > oo before the continuum limit is taken |]l8| . 
Thus, DWF fermions have the full axial symme- 
try of the continuum at any lattice spacing. How- 
ever, this does not rule out the possibility that the 
doublers reappear at strong coupling in order to 
trivially satisfy Eq. (||). For the singlet case, the 
last term gives rise to the usual axial anomaly |]l9|. 

3. NUMERICAL RESULTS 

The first numerical works using DWF for vec- 
tor gauge theories were studies of the two dimen- 



sional vector Schwinger modcl||2^,|2^. Both indi- 
cated DWF were practical and useful for numer- 
ical work. An in-depth study of DWF using the 
two dimensional vector Schwinger can be found 
in 0. 

The first simulations of QCD using DWF are 
found in Refs.||lO|-[l^. There it was shown that 
even modest numbers of sites in the extra dimen- 
sion were sufficient to maintain the desirable chi- 
ral properties of DWF. Within the last year a 
significant interest in using DWF for QCD has 
developed, leading to many new results that were 
presented at this conference. These are described 
below. 

3.1. PCAC Ward identity, zero modes, 
topology, and finite Ng effects 

The anomalous term in Eq. (^) for the PCAC 
Ward identity (O = J5) can be used to estimate 
the effect of finite Ng in simulations. Taking the 
ratio of Eq. (^) with the pseudoscalar correlator 
leads to 2m+2 {J5q{x)J5{y)) / {J^{x)J^{y)) on the 
r.h.s., so, as long as the second term is small com- 
pared to the explicit quark mass m, one expects 
good chiral behavior. 

Results for the relative anomalous contribution 
are shown in Fig. (|l|). For M — 1.7 all the points 
are at least an order of magnitude less than the 
value of m(0.05) used in the simulations, except 
the point at Ng = 10 and (3 = 5.85. As shown 
below, these simulations display good chiral be- 
havior whereas the pion mass, for example, at 
(3 = 5.85 and Ng = 10 has a small but non-zero 
value when extrapolated to m = 0. Note also that 
the effects decrease significantly as the coupling 
is weakened while the rate appears to increase. 

The data in Fig. (|l|) is plotted on a semi-log 
plot, so the four points at {f3 — 5.85, M = 1.7) 
should lie on a straight line if the suppression 
of chiral symmetry breaking effects is given by 
a simple exponential. However, there is notice- 
able curvature (but, note that the sample of con- 
figurations is small) . This curvature may indicate 
power law suppression or exponential suppression 
with an Ng dependent rate. The latter was ob- 
served in the two dimensional vector Schwinger 
model [0 and was attributed to gauge fields hav- 
ing trivial versus nontrivial topology. 
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Figure 1 . The relative anomalous contribution to 
the PCAC Ward identity. 

The former case may arise from translationaUy 
invariant modes in the 5th dimension which are 
expected to be the dominant anomalous contri- 
butions to the Ward identities Such modes 
may arise when the 4d Wilson-Dirac operator 
{pw{—M)) supports an exact zero mode(here 
—M refers to supercritical values of the hopping 
parameter) . This is because the spectrum of the 
Hamiltonian Hw{—M) = j^w{—M) defines a 
transfer matrix for propagation in the 5th dimen- 
sion§, = eM-Hw)^'- When Hwi^M) 
has a zero eigenvalue, T has a unit eigenvalue, and 
there is no suppression in the extra dimension. 
There is a gap in the spectrum below the M cor- 
responding to Kc, which vanishes at Kc- Naively, 
the gap is expected to reopen above this value 
of M|2|. The spectrum of Hw{-M) has been 
calculated recently |2^,|2^ on quenched configura- 
tions corresponding to those in Fig. (|l|). It was 
found that Hw{—M) contains many zeroes, or 
level crossings, arising from lattice instanton-like 
artifacts, which could close the gap for all values 
of M suitable for DWF simulations, a disaster. 

The data for DWF, on the other hand, are not 



consistent with this scenario. If there were no 
suppression, light modes would not be bound to 
the walls and the ratio of pseudoscalar densities 
in Fig. (|l|) would be 0(1). However, we have al- 
ready seen that this ratio is small for M — 1.7 
which is far from the Kc region. For M — 1 
which is near this region at /3 = 5.85, the ra- 
tio increases dramatically, signalling the absence 
of light modes. This is consistent with the gap 
closing at Kc, as it should. It may be somewhat 
confusing that the existence of light modes for 
DWF rules them out for Wilson fermions, and 
vice- versa. 

TranslationaUy invariant modes do exist in 
DWF simulations, however, so it is important to 
understand their impact on the 4d physics. In 
Fig. (||), I show ip-tp as a function of the extra 
5th coordinate. The plot is from the Columbia 
group and is for a single gauge configura- 
tion at /3 = 5.85. The physical 4d value is at 
s — Ns — I. M — 1.45, which corresponds to an 
exact zero of Hw- One sees clearly the invariant 
mode, but its contribution to the physical value of 
'0V' is only about one percent. Note, zero modes 
of Hw{—M) are not zeroes of DWF. In fact, these 
modes may be quite heavy; since mixing between 
the modes on opposite walls is maximal in this 
case, the induced mass may be large. 

The Columbia group has also calculated (00) 
on an ensemble of 200 gauge fields at quenched 
f3 — 5.85 In Fig. (||) a divergence as m ^ 
is clearly evident. The divergence arises in the 
ensemble average from zero modes that are un- 
suppressed by the (missing) fermion determinant. 
This has not been seen in lattice simulations un- 
til now, presumably due to lattice artifacts. The 
coefficient of the l/m term should decrease by 
1/VV. The results in Fig. (|) are for V = 8^. 
The l/m term decreases by a factor of six on a 
F = 16^ lattice. The Columbia group is currently 
running on a still larger lattice to resolve this vol- 
ume dependence. 

Finally, the Argonne group has been calculat- 
ing the low lying spectrum of Rj^Ddwf in order 
to study the behavior of fiavor singlet correlation 
functions and the restoration of the U{1) axial 
symmetry above the chiral symmetry phase tran- 
sitionpffl. A sample of their results is shown in 
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8 x32, |3=5.85, quenched, 200 configs 
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Figure 2. ^pip as a function of the extra coordinate 
s on a single quenched configuration at /3 = 5.85 
from the Columbia group The value of M 
corresponds to a zero mode of Hw{—M). The 
physical value is on the boundary, s — Ns — 1. 



Fig. (Q) for quenched /3 = 6.2. Similar results 
hold at /3 = 6.0. Note that TV^ < 10 in the figure. 
The lowest eigenvalue is clearly exponential in Ns 
while the next lowest is iV^ independent. For 
configurations with a net number of instantons, 
they find the same number of exponentially small 
eigenvalues, in accord with the Atiyah-Singer In- 
dex theorem. They explicitly calculate the topo- 
logical charge Qtop = ™/ d'^xq^'j^qx (Fig. (||)) 
and find that it is integer valued within small er- 
rors, an impressive result considering the values 
of Ng used. The deviation of Qtop near the ori- 
gin is a finite effect. Strictly speaking, the 
Index theorem for DWF is only exact in the limit 
Ns — > oo. Away from the continuum DWF (and 
the overlap) do not have a unique index; it de- 

for example). 



pends on M(see Ref. |27 



Thus, I conclude that the gap in Hw{—M) is 
not actually closed for DWF simulations at val- 
ues of the coupling used in current simulations. 
Instead, DWF simply support a number of zero 
modes corresponding to the net number of cross- 
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Figure 3. (■(/'V') from the Columbia group p5| as a 
function of the quark mass. The presence of un- 
suppressed zero modes gives rise to a divergence 
as m ^ 0. 



ings in the spectrum of Hw below a chosen value 
of Af , as they should. In the ensemble average 
and y — > cxD, as long as the density of exact zero 
modes of Hw is vanishingly small for a particular 
M, then M is suitable for DWF simulations. 

There is a trend for the suppression to weaken 
as (3 decreases(see Fig. (|l|)). In Ref. it was 
found that Ns = 18 for (3 = 5.85 was necessary 
for the pion mass to extrapolate to zero within 
statistical errors compared to Ns — 10 at 6.0. 
At 5.7 the Columbia group found that even for 
Ns = 48, the pion mass has a significant in- 
tercept at m — OH^. However, it is not clear 
that this is entirely a finite Ng effect since results 
from Ns — 24 are roughly the same. The inter- 
cept may be related to unsuppressed zero modes 
from quenching. Our experience has been that at 
strong couplings there is a minimum value of Ns 
as well, and tuning M does not have a large effect. 
The above behavior may be related to an Aoki 
phase of Wilson fermions at strong coupling p8|. 

3.2. Light quark masses and spectrum 

The quenched QCD pseudoscalar spectrum us- 
ing DWF was first studied in Refs.|p^|-[l^, along 
with rough estimates for the vector channel. 
These measurements have also been used to calcu- 
late the strange quark mass which was presented 
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operator vs. Ng on a configuration with Qtop = f • 
The plot is from the Argonne group p^. 
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Figure 5. The topological charge Qtop from the 
Argonne group, Rcf. [p6| , on a single configura- 
tion with one anti-instanton. 



at this conference |29|]. The Columbia collabora- 
tion presented quenched results for the light pseu- 
doscalar, vector, and nucleon channels at /? = 5.7, 
as well as the chiral condensate. 

In Fig. the pion mass squared is plotted as 
a function of m for 1.5 < M < 2.1 at f3 — 6.0. 
Each curve linearly extrapolates to zero, in accord 
with lowest order chiral perturbation theory. Fig. 

shows the same data versus M. has the 
characteristic quadratic shape expected from free 
field theory but with a maximum that is shifted 
from its tree level value. From the figure the shift 
is approximately 0.8, in good agreement with the 
estimates of Mc above. The dashed lines corre- 
spond to fits to Eq. (H) with Mc fixed to the 
value corresponding to Kc for Wilson fcrmions(see 
Ref. H for details). 

The vector and nucleon masses from the 
Columbia group are shown in Fig. (H). The 
quenched coupling is 5.7 and the lattice volume 
is 8'^ X 32 X 48, so these correspond to roughly 
the same physical volume as the above(in fact, 
the rho mass scales within errors with the result 
in Ref. [0). From their fits, mnuc/nT-p ~ 1-4 in 
the limit m = compared to roughly 1.5 calcu- 



lated using Kogut-Susskind quarks at the same 
physical volume ||30[|. As mentioned above, the 
chiral symmetry breaking effects due to finite Ng 
are larger at stronger coupling, so the Columbia 
group has used TV^ = 48. Even so, the pion mass 
is not zero in the limit to = 0, but corresponds to 
roughly 230 MeV. 

The Columbia group has calculated (i/'V') for 
several values of M at quenched (3 = 5.85 (Fig. 
(^) which yields a clear demonstration of the fla- 
vor structure of DWFQ. Recall that in the 
free theory we had one flavor for < M < 2, 
four flavors for 2 < M < 4, and so on. A simi- 
lar pattern emerges nonperturbatively, once again 
shifted to larger M. The light quark (■(/'V') first 
becomes non-zero at Mc, flattens out in the region 
1.65 ^ M ^ 2.15, rises rapidly again, and turns 
over at about 3.4. The first region corresponds to 
one flavor and the second, four. The ratio of {ijjilj) 
in the two regions is numerically close to four. In 
the free theory the regions are quadratic polyno- 
mials with zeroes at 0, 2, 4, ... , 10. Note that the 
results become Ng independent at Ng = 24. 

A first estimate of the (quenched) strange 
quark mass was presented at the meeting p^. 
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Figure 6. The pion mass squared as a function of 
the quark mass. Ng = 14, /3 = 6.0. 



The ^1 — 2GeV result is compared to conven- 
tional determinations in Fig. (^0|), where one sees 
rough agreement. The DWF result also appears 
to scale nicely with the lattice spacing, albeit 
within rather large statistical errors. A weighted 
average gives nis = 82(15)MeV. 

To get rus , several steps are needed. First m is 
tuned in the usual way to yield the kaon mass. To 
match the DWF result with a continuum scheme, 
the one loop self-energy calculation of Aoki and 
Taniguchi was extended to the massive case, m 7^ 
0. One must determine a value for M = M — Mc 
to stick into the one loop analog of Eq. (^) . In the 
above it was determined from Kc- This method 
is not exact, even though it is nonperturbative. 
That is because of the symmetry of the action 
under the change M ^ 10 — M; M = 5 is a fixed 
point so M is not shifted uniformly. While a more 
accurate determination of M is desirable, simply 
using Kc only introduces a small error in the final 
value of M(2-M)Q. Ultimately, the best solution 
is to calculate the quark mass renormalization Zm 
nonperturbatively (using the method of |^l|, for 

^ I thank Y. Shamir for pointing this out. 



Figure 7. The pion mass squared as a function of 
the domain wall height. The plot is from Ref. [ p9[ . 
The dashed lines correspond to fits to Eq. ^) 
using M and the nonperturbative value of Mc- 



instance), where this extra factor is taken into 
account automatically. 

3.3. Weak interaction matrix elements 

The study of weak interaction matrix elements 
was one of the original motivations for using DWF 
for QCD simulations [ p^|jTl| ] . This is because the 
matrix elements of many important weak opera- 
tors between pseudoscalar states vanish linearly 
with the quark mass, which is a direct conse- 
quence of chiral symmetry p^. Ordinary Wil- 
son quarks explicitly break chiral symmetry, so 
the naive lattice transcription of these operators 
does not vanish in the chiral limit. This is a well 
known, long standing problem. In principle, the 
solution is to fine tune the continuum-like op- 
erators with others that transform under differ- 
ent chiral representations to restore the correct 
behavior in the chiral limit. Alternatively, one 
may use Kogut-Susskind quarks which do exhibit 
the correct behavior since they maintain an exact 
abelian axial symmetry on the lattice. This rem- 
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Figure 8. The vector and nucleon masses as a 
function of m. The plot is from Ref. 



nant of the continuum is enough to ensure the 
proper behavior. However, the price of Kogut- 
Susskind quarks is the breaking of the continuum 
flavor symmetries which induces significant lat- 
tice artifacts. As we shall see, DWF seem to have 
significantly smaller discretization errors than the 
Kogut-Susskind variety. 

The vanishing of weak matrix elements in the 
chiral limit is a stringent test of the chiral symme- 
try of DWF since even with operator (un) mixing 
this is a difficult achievement (see Refs.[ ^5||33| for 
recent results with nonperturbative and pertur- 
bative mixing, respectively). Fig. ( pT] ) shows 
the matrix element of the left-left weak operator 
{Oll = (s7p(l - 75)d)^) that describes K° ~ k° 
mixing. For all three values of /3, 5.85, 6.0, and 
6.3, it vanishes linearly with m, as it should ac- 
cording to chiral perturbation theory [Q. 

The matrix element of Oll is conventionally 
normalized to its value in the vacuum satura- 
tion approximation which yields the B parame- 
ter. The kaon B parameter, Bk, is an impor- 
tant phenomenological parameter used to extract 
part of the CKM matrix from experimental mea- 
surements, which in turn constrains the Standard 
Model. The quenched DWF result is shown for 
several lattice spacings in Fig. The impres- 

sive Kogut-Susskind result of Ref. is shown 
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Figure 9. (V'V') as a function of M from Ref. p4| ]. 



for comparison. Since Oll is generated by the 
operator product expansion, it is scale dependent, 
and therefore the result at each coupling must be 
run to a common scale before the continuum limit 
is taken. In the continuum, the vacuum satu- 
rated value of {K^\Oll\K'^) is proportional to the 
square of the axial current matrix element which 
does not get renormalized since the axial current 
is partially conserved. As explained above, DWF 
have a partially conserved axial current, but this 
is not the current that is needed for the vacuum 
saturation of Oll on the lattice since we have 
chosen to construct Oll from the quark operators 
defined in Eq. (^). Therefore, the denominator 
in Bk also receives a (finite) multiplicative renor- 
malization. Since these renormalizations have not 
yet been done for DWF, the most that can be said 
is that the result for Bk is consistent with the 
Kogut-Susskind result. Nevertheless, the indica- 
tion from Fig. ( [l^ ) is that discretization errors 
are significantly smaller in the DWF case. Need- 
less to say, it is a high priority to reduce the sta- 
tistical errors and determine the renormalization 
factors in order to precisely determine the lattice 
spacing dependence. It is important to point out 
that the lattice volumes of the last two Kogut- 
Susskind points in Fig. jl^ ) are 56'^ x 96 and 
48'^ X 96 and correspond to physical volumes of 
2.3fm^ and 2.8fm^, respectively. If the improved 
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Figure 10. The strange quark mass from Ref. 
after matching to the MS scheme. 



Figure 11. The matrix element of Oll between 
degenerate pseudoscalar states. It vanishes in the 
chiral hmit, a stringent test for the chiral symme- 
try properties of DWF. 



scahng of DWF holds up, then the added cost of 
the extra dimension will be significantly offset by 
the smaller Ad lattice volumes allowed by larger 
lattice spacings. 

At f3 — 6.0, there is also a value of Bk on a 
24'^ X 40 lattice which indicates that finite volume 
effects are not large. 

Encouraged by the results of our Bk calcula- 
tion, we have begun calculations of the As = 1 
effective weak Hamiltonian operators that govern 
K TTTT decays and are necessary for e'/e and the 
AJ = 1/2 rule, for example [p^. Once again, in 
the continuum, chiral symmetry restricts the be- 
havior of these operators and is therefore impor- 
tant for the lattice calculations as well. In fact, 
the good chiral properties of DWF allow the use 
of the method of Bernard, et al. to calculate 
unphysical K n amplitudes and relate them, 
using chiral perturbation theory, to the physi- 
cal iiT — > TTTT amplitudes. The unphysical three 
point functions are easier to calculate on the lat- 
tice. Kilcup and Pekurovsky are using the same 
method with Kogut-Susskind quarks to calculate 
these amplitudes jSTj . There it is already clear 



that perturbation theory fails for the renormal- 
ization of these operators, and a nonperturbative 
method is needed. On the other hand, since the 
one-loop renormalization of the DWF quark mass 
is reasonable, it is still possible to renormalize the 
DWF operators perturbatively. 

Our preliminary results on an ensemble of eigh- 
teen gauge configurations at (3 = 6.0 are shown 
in Fig. The lattice volume is 24^ x 40 x 10. 

As is well known, these amplitudes are more dif- 
ficult to calculate than Oll because of the no- 
torious "eye" contractions which occur since a 
quark and an antiquark with the same flavor in 
the operator can annihilate each other. In order 
to efficiently average the operator over all spa- 
tial sites on a fixed time slice, one uses a random 
source(s) on the time slice. The results shown in 
Fig. (^_3|) are for a single random source in the 
middle of the lattice. It is therefore encouraging 
to see a signal at all. Oi, . . . ,Oq which transform 
under (8,1) representations of SU{3)l x SU{3)ji 
should vanish linearly in the chiral limit. O5 and 
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Figure 12. The kaon B parameter. The Kogut- 
Susskind result is from Ref. |3^. DWF indicate 
improved scahng in this case. 



Figure 13. K t: matrix elements of the As — 1 
effective weak Hamiltonian. 



Oq do not, though they miss the origin by less 
than two standard deviations. This could be due 
to the small statistics of our ensemble, or that 
Ns = 10 is not large enough. Also note that 
there are only two quark masses, and the ran- 
dom source was different for each. O7 § trans- 
form under (8,8) representations so do not vanish 
in the chiral limit ||3^ . 05,...,8 are expected to be 
larger than Oi_...^4 since they are left-right oper- 
ators which is the case here, and the color mixed 
operators are roughly three times their unmixed 
counterparts, as expected from the vacuum satu- 
ration approximation. 

In order to carry the above through to the phys- 
ical K —f TTTT decays, an important vacuum sub- 
traction, {K\Oi\0) must still be madep^. 

3.4. Dynamical simulations and non-zero 
temperature 

The Columbia group has begun the first simu- 
lations using dynamical DWF to study the two 
flavor QCD chiral symmetry restoration phase 
transition at non-zero temperature. These sim- 



ulations are obviously quite demanding. Impres- 
sive preliminary results were presented at this 
conference |l^. They, as well as the Argonne 
group also presented quenched and semi-classical 
results related to the anomalous t/(l) axial sym- 
metry [^,^ 26 1 . A detailed description of the 
Columbia group's semi-classical results can be 
found in Ref. |^. 

Using DWF, the two flavor chiral phase tran- 
sition in QCD appears qualitatively similar to 
the Kogut-Susskind result. In Fig. ( p!^ ) there 
is a smooth crossover in both the Polyakov loop 
and (ipip) starting at /3 « 5.2 {N^ = 4). In the 
limit m ^ 0, below the transition (V'V') is sponta- 
neously broken while above the transition chiral 
symmetry is restored (Fig. (15)). The very small 
but statistically significant non-zero intercept in 
Fig. ( |l5| ) is presumably due to finite Ns effects. 
From Fig. ( p^ ) ordinary finite volume effects are 
small as well. While these results are very en- 
couraging, the quark masses are still too heavy 
to discuss details like the order of the transition 
and critical exponents. It should be stressed that 
for the first time the correct number of light pi- 
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ons is being simulated in dynamical studies of the 
Uf — 2 QCD chiral phase transition, which may 
be crucial to the dynamics. 
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Figure 14. (ipip) and the Polyakov loop for Uf ~ 2 
from the Columbia collaboration [l^ . Both ex- 
hibit a rapid crossover near /3 = 5.4. 
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Figure 15. ('0V') as a function of m for /3 = 
5.20(upper curve) and 5.45(lower curves). Chiral 
symmetry is restored above the phase transition. 
The two lower curves correspond to 8'^ and 16^ 
lattices. Uf — 2. The plot is from the Columbia 
group Ref. M. 



The pseudo-critical coupling is a function of 
M(Fig. (p^). There are two reasons for this: 
the quark mass depends on M (Eq. (H)), and 
the number of flavors is a function of M. Below 
Mc there are no light states bound to the domain 
wall. At Mc the domain wall states just begin 
to form so f3c starts at its quenched value. As 
M increases, a full range of momenta are avail- 
able to the states on the walls and /3c decreases 
to its = 2 value(remember that the fermion 
determinant is D^D). As M increases still fur- 
ther, four light quarks eventually appear on the 
domain walls so /3c decreases to its = 8 value. 

Both the Argonne and Columbia groups have 
studied the restoration of the U{1)a symmetry 
above the chiral phase transition. As mentioned 
earlier, the Argonne group has measured the 
topological charge(see Fig. (^)) on quenched con- 
figurations, but have not yet published their re- 
sults for flavor singlet meson masses, correlators, 
or susceptibilities. The Columbia group has mea- 



sured Uf = 2 screening masses for the pion and 
the delta which should be degenerate if the U{1)a 
symmetry is restored. The mass difference some- 
what above the chiral transition is shown in Fig. 
(p^). There is a small but significant value at 
m = 0. Note the data appears quadratic in m, as 
expected, due to the restoration of chiral symme- 
try. There are no known lattice artifacts present 
to obscure this picture as is the case for Kogut- 
Susskind quarks. 

In general, dynamical simulations require a 
Pauli-Villars(PV) subtraction. The action for the 
Ns — I heavy states induced by the extra di- 
mension dominates the total action in the limit 
Ns ool^. The form of the subtraction is not 
unique; a second order operator is suggested in 
Ref. 1^ while in Ref. 0Jl^ a first order operator 
is used. The latter version is vulnerable to the 
"exceptional" configurations discussed in section 
3.1. That is, if Dw{—M) supports an exact zero 
mode, so does the first order operator, which has 
periodic instead of anti-periodic boundary condi- 
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Figure 16. The Polyakov loop from Ref. [|6[ for 
several values of M. nj — 2. The pseudo-critical 
couphng depends on the choice of M. 

tions in the extra dimension. This zero mode is 
not canceled by a corresponding massless mode in 
the DWF operator. On the other hand, the au- 
thors of Ref. ||l6| have seen no evidence for these 
modes in their simulations. A method for sup- 
pressing such configurations has been suggested 
in Ref. |40|, which is described in the next sec- 
tion. 

4. TOPICS FOR FUTURE STUDY 

Beyond using DWF in their present form to do 
useful phenomenology, future studies will focus on 
improvements that reduce the size of the extra 
dimension and/or increase the rate of exponen- 
tial suppression. Several proposals have already 
been made ||4l| , ^ , ^ . The first proposal uses a 
hyper-cubic action instead of a Wilson like action. 
The last two include an explicit lattice spacing in 
the extra dimension which renders the DWF ac- 
tion positive definite for any gauge coupling. In 
Ref. |4^ it is argued that this restricts eigenval- 
ues of the transfer matrix on the unit circle to the 
real axis, reducing slowly decaying chiral symme- 
try violations. There, a modified PV action is 



0.4 



0.3 



0.2 



E 

E 0.1 



-0.1 



16^x4, p=5.45, M=1.9, Nj=16 




0.05 



0.1 

m 



imercept = 2.1(1.3)x10 

0.15 0.2 



Figure 17. The difference in screening masses 
above the chiral symmetry restoration phase tran- 
sition from Ref. |16|. Uf — 2. A difference of zero 
at TO = indicates restoration of the U{1)a sym- 
metry. 

considered which suppresses the unwanted zero 
modes of D\y{—M); an extra term involving the 
4d links is added to the PV action to do the job. 
A side-effect is the small renormalization of the 
coupling constant. 

An exciting development is the related deriva- 
tion of the overlap-Dirac operator by Neu- 
berger[0. This lattice discretization, which cor- 
responds to a single massless Dirac particle, is 
derived from the overlap, but retains no refer- 
ence to an extra flavor space or dimension. It is 
completely 4d and has been shown explicitly to 
satisfy the Ginsparg- Wilson relation|^^. Initial 
studies indicate the method may be feasible nu- 
merically Q-Q, though its application to lattice 
QCD has not yet been demonstrated. 

Finally, a new approach to lattice gauge the- 
ory is the quantum link model. DWF arise natu- 
rally in quantum link models since they are con- 
structed in five dimensions . 

5. CONCLUSIONS 

I have reviewed the present status of DWF 
for vector gauge theories, in particular QCD. We 
have seen that DWF with a modest number of 
sites in the extra dimension maintain their good 



13 



chiral properties, which are exact in the limit 
A'^s — > oo. The results show promising behavior in 
studies where chiral symmetry is important such 
as weak matrix element phenomenology, dynam- 
ical simulations of the chiral symmetry restora- 
tion phase transition, and the relation between 
topology and fermionic zero modes. Indications 
of improved scaling with DWF need to be verified 
since this may significantly offset the cost of the 
added dimension. 

I acknowledge many helpful discussions with 
M. Creutz, R. Mawhinney, Y. Shamir, A. Soni, 
and M. Wingate. This work was supported by the 
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